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1 Introduction 



Measurements of high energy exclusive production of J/ip meson in photon-proton colhsions 
serve as an important testing ground to quantify the "hard" physics of pQCD in the hmit of 
a small scaling variable x. A process is considered to be "hard" if a large momentum scale is 
involved, so that the leading contribution of the hadronic fluctuation of the photon is a small 
perturbative qq pair. This large momentum scale may be a heavy quark mass, large virtuality 
of the photon or large four momentum transfer squared at the proton vertex (— t). For the 
diffractive photoproduction of J / ip this hard scale is provided by the charm quark mass. In our 
recent letter |I| we presented a screened model which yielded impressive agreement with the 
rich experimental data on high energy J/ ip photoproduction. The model includes contributions 
of the off diagonal (skewed) gluon distributions 0, as well as the real part of the production 
amplitude. In addition, we addressed the issue of Fermi motion of the heavy quark within the 
quarkonium system which provides further suppression of the final cross section. We argued 
that on one hand this effect is sensitive to the charm quark mass, while on the other hand, in the 
relevant kinematical region of the experimental data, it has a rather mild energy dependence. 
An obvious conclusion of the above argument is that with the current theoretical uncertainties 
regarding the Fermi motion suppression, the benefits resulting from a phenomenological model 
would be rather small. The correction due to Fermi motion has been estimated as an overall 
suppression factor of approximately 30%. 

In the present work, we present an analysis for the t dependence of the J/ ip differential cross 
section which is compatible with our recent investigation of the integrated cross section. To 
compare with the experimental data we need to know B - the J/ ip forward slope of the differ- 
ential cross section. In [|l| we assumed a weak dependence of B on x, which was calculated by 
averaging the square of the impact parameter b, using our screening correction (SC) formalism 
[If. This input assumption is not sufficient for a detailed t dependence study of J/ip differential 
cross sections as the presently available preliminary experimental data on i? Q show a steeper 
energy dependence than our prediction. The goals of the present improved version of our model 
are, thus, to reproduce the value as well as the energy dependence of B while maintaining the 
quality of the results achieved in 

In our previous study, a 6-space Gaussian was assumed for S{b), the gluon ladder-proton 
non perturbative form factor |]T|, ^, Q . After a careful check, it appears that the integrated cross 
sections and nucleon structure function [|, |^ are not sensitive to the choice of S{b). The 
impact parameter dependence plays an important role, when calculating screening corrections 
to heavy nuclei structure functions, where a Wood-Saxon like profile function has been used 
0, or for the J/ip differential photoproduction cross sections [|lOl which is the main goal of the 
present study. As we shall demonstrate, the experimental data |^ for the forward differential 
cross section slope can be well reproduced with the observation that the energy dependence 
of B is sensitive to the choice of S{b). It appears that the Fourier transform of the dipole 



1 



electromagnetic form factor, is the preferred profile for estimating the screening corrections to 
B. 



The momentum transfer t at the proton vertex controls the size of the interacting system, 
which emits the gluon ladder. The t dependence of hard processes can, thus, be separated into 
three kinematical zones. The low |t| region, where the transverse size of the interacting system is 
of the order of l/Qo; where Ql is the input scale for the DGLAP|]TT| evolution {Ql ~ 1 GeV^). 
This region, in which < \t\ < Q^, has been explored experimentally, mainly via elastic 
processes |T^. The hard scale in this region is determined by the c-quark mass ~ 
The second region is the intermediate region Qo < |^| < M^. We suggest that the physics of this 
region is similar to the physics of the small \t\ region, i.e. the process is still dominated by large 
logarithms due to momentum integration over loops in the ordered gluon ladder. However, 
the scale of the initial conditions of the DGLAP evolution is determined by t, as it becomes 
of comparable size to Q^. More specifically, the process is dominated by \og{Mj^^/t). In the 
intermediate region, the vector meson is produced mostly quasi elastically, while the proton 
dissociates into a diffractive mass. Note that the diffractive t slope is considerably smaller than 
the elastic one. In the present study we compare our calculations in this region with preliminary 
experimental data which were read off figures given in Refs. |]13|, As |t| increases, the hard 
process behavior is dominated by the value of the \t\ dependence is inverted [jl5| and the 



momenta on the gluon ladder are no longer strongly ordered. Consequently, predictions for the 
cross section cannot be described using conventional DGLAP evolution. This third region of 
large \t\ has not been explored experimentally. 



These three regions of |t| were first discussed in |]T^ by Forshaw and Ryskin and by Bartels, 
Forshaw, Lotter and Wuesthoff in the framework of the BFKL equation [0, where it was 



stated that each of the above regions corresponds to different physics. In a later publication 
ir^ , a BFKL calculation of vector meson production, and in particular of J/ip production, was 



presented and a good description of experimental data was obtained. 

In this paper, we suggest a description of the second kinematic region of intermediate \t\ 
{Qo < 1^1 ^ using the DGLAP [|ll|] approach, this yields simple and transparent formulae 
which have explicit matching to the region of small \t\. 

Our paper is organized as follows: in section ^ we calculate the forward slope for different 
profiles and compare with the experimental data. In sections ^ and ^ we calculate the differen- 
tial cross sections for exclusive J/ip production both in the elastic and the inelastic channels, 
respectively. Our conclusions are given in section 5. 



2 Forward Slope for Differential Cross Section 

In this section we calculate the 5-slope for J/ip production, including screening corrections for 



up to one extra gluon emission. Our formalism for calculating the screening corrections |18 
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is based on the iteration of the non-hnear evolution equation |[T9[ for the imaginary part of 
the elastic amplitude for a dipole to scatter off a hadron target, which is valid in the whole 
kinematic region, including the low x region EO . 



The interaction of a dipole with the target is realized through an exchange of gluons. In 
the target rest frame, one can use Gribov's factorization in the context of a dipole-nucleon 
interaction: a dipole emits a soft gluon which develops into a ladder by successive emissions 
of small X gluons. Subsequently, the gluon ladder interacts with the target. At high energies, 
the transverse size of an interacting dipole does not change during the QCD interaction, and 
therefore dipoles are good degrees of freedom. In the limit of large number of colors, a leading 
logarithmic soft gluon wave function is equivalent to a dipole wave function, and thus, the 
interaction between a dipole and a target can be treated as a process which is subsequent to a 
transition of a dipole into two dipoles. Consider a decay of a dipole with transverse size r± into 
two dipoles as shown in Fig. |I|. The probability for this decay is given by the square of the wave 
function of the dipole. At large Nc one can view this decay as an emission of a zero transverse 
size gluon. In this framework, the cross section for the virtual photon-nucleon interaction is 
written as 

(^totil*p) = J d^r^ j dz ■^^'{Q'^-^r^^z) ^ adipoie{r±, z) (1) 

where the wave functions of the virtual photon are well known, in particular the transverse 
wave function is given by the Bessel function Ki. In our model, the dipole-nucleon cross 
section consists of two contributions: 



The first, denoted a'^{r±), is the percolation of a qq pair through the target, without taking 
into account the extra gluon, and it is calculated using a multi gluon ladder exchange. 
This contribution can be expressed as a product of the dipole wave function and the cross 
section for the dipole to scatter off the nucleon. In impact parameter space, a"^ can be 



written in the form ||T8|, 21 



a'{r^) oc J d%^(fri_Ki{Qr^) (l - e~5f^«) , (2) 

where Q = {Q'^ + M^) /4, x = {Q^ + M^) /W^ and the opacity VL^ = Qg^b, r_L, x) is defined 
as 

= Vlq(b,r^,x) = y^l"s (^-^^ xG (^x, -^^ S{b±). (3) 

The b± dependence of Qq is expressed through the profile function S{b±), which we shall 
discuss later. Note that the argument of the Bessel function should be ar±, where = 
z{l — z)Q'^ + however, for heavy quarks, we use the approximation z = \. 

For the second contribution, denoted (3"^(rx) below, we consider an additional gluon in the 
hadronic state which percolates through the target. We assume that the parent dipole 
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Figure 1: A decay of one dipole with transverse size r± into two dipoles. In the case where 
two of the produced quarks are with smah transverse separation, they can be view 
as a zero transverse size gluon, hence the final sate is a qqg state. 



whose transverse size is r±, is larger than the produced dipole whose transverse size is p± 
(see Fig. |lD. The second contribution can be written in the form [l^, |2T[| : 



dx 



d^pL 



(1 



e 



(4) 



with ng{b,r±,x) = r^, x). 



The forward slope B is given by [p2 

B 



1 , _ JbldblSjb^) 
2^^' 2jdb\S{b^) 



(5) 



where the factor of one half is introduced since B is measured for cross sections rather than 
amplitudes. To calculate the forward slope -B, we need to average over the impact parent 6^, 
using the sum of @j and @ as weights, 



B 



Jd%^d\^blKi{Qr^) 






2Jd%^d^r^Ki{Qr^) 


(l-e-^^') 





(6) 



Being the root mean square average in b± space, B is sensitive to the b± dependence of the 
opacities Qq and Qg, which depend on S{bj_). 

The profile functions S{b±) are connected to the form factors via Fourier transforms. We 
have investigated three different form factors listed below, and compared the results with the 
relevant experimental data on the J/ forward slope . 
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1. An exponential form factor, 

Fe.,it) = ei^'*, (7) 
which transforms into a Gaussian in impact parameter space. 



2. An electromagnetic form factor 



Fdipole{t) ~ Q _ lp2^^2' 



its Fourier transform in b± space is 



3. For completeness, we have also investigated the form factor, suggested in p3 



4ml- 2-79^ 1 
where rup is the proton's mass. The Fourier transforms of F^l is given by. 
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S{b±) = — ai (Ko{Vo:71b^) - Ko(2mp6^) ) + a2bKi{V0Jlb^) , (11) 

where the values of ai ~ 0.8 and 02 ~ 0.33 were calculated from the numerical coefficients 
in fllOD. 



The various coefficients are determined from the normalization condition / d?b±S{b^_) = 1. 

In our calculations we have substituted each of the above mentioned profiles in (^), and 
calculated B for the three different PDF parameterizations [^, ^ for the gluon distri- 
bution which appears in (j^). Our results are shown in Fig. H. While the energy dependence 
obtained using Fdipoie and Fdl is in agreement with the ZEUS data, the Gaussian distribution 
predicts an increase with energy which is too mild to reproduce the measured slope. Regarding 
normalization, the data is well reproduced by (P) with the choice = 10 GeV~^. However, 
some adjustments are needed regarding the numerical parameters of (]TT|). The choice of our 
single parameter, is comparable with our earlier estimations (i?^ = 8.5 GeV~'^), where we used 
a Gaussian distribution for S{b±), and we maintain the quality of our results presented in 
Ref. m. 
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Figure 2: The forward slope for differential cross section of J /ip production, data and pre- 
dictions, corresponding to different profile functions S{b±). Predictions have been 
calculated for three different PDFs. 
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3 Differential Cross Section at Low It 



In this section we calculate da/dt as a function of t for the region of \t\ < Qq, where Qq is 
the scale at which the input for the evolution is computed. Our expressions are based on a 
simplified model of two Pomeron exchange, which we find to be satisfactory, given the sizeable 
experimental errors, although our screening correction formalism is capable of deriving a full 
series calculation. It should be appreciated that what follows in this section is a very simplified 
picture of screening corrections to the lowest order and it should not be used for calculating 
other physical quantities, e. 5^. integrated cross sections, on which more accurate data exist. 

The procedure for calculating the forward differential cross section for photo production of 
a heavy vector meson in the color dipole approximation is straightforward P^. We follow |I| 



and write the differential cross section at t = 0, in the leading logarithmic approximation of 
pQCD, including a contribution from the real part of the production amplitude and the 
skewed (off diagonal) gluon distribution, 

MJP^VP) _ 16.'r„ ^2-«r(A + i)f /l + ta„^(^)| ai(Q>G^=^,Q\ (12) 




dt SctemM^ \ ^^F^fA + 4) ) \ '2 



where A = (91ogxG'/(91og(l/a;). In the derivation of (12), a simple static non relativistic 
estimate of the vector meson wave function was used. As stated we consider the correction 
due to Fermi motion, which changes the numerical value obtained from the above equation, by 
an overall (constant) suppression factor. 



In [|T| we took into account further suppression of (|T2D due to screening. These are facilitated 
by damping factors |^, ^ which are the ratios of the screened to non screened observables. 
Specifically, the DLLA contributions to the damping factors, are due to screening in the quark 
sector (denoted Dg) and and in the gluon sector (denoted Dg). The DLLA screening correction 
calculation, was derived in pOf using a Gaussian distribution for the profile S{b±) [see Eq. 



For completeness, we briefly described the basic formulae. 

The damping factor due to the screening in the quark sector i.e. the percolation of the cc 
through the target, is: 



where k„ is given by 




Di = -Ei(-)e'=« , (13) 



= ^xG^^"'^^(x,Q^). (14) 

Our expression for the damping in the gluon sector, is the square of the gluon damping defined 
in 1,0], 

xG'^ix,Q') = D,(x,g2)xG'^^^^^(x,g2), (15) 
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Figure 3: The first two contributions to the amplitude. 



where 

71'' Jx X' Jo J ^ ' 



and Kg = jKq. The overall damping factor for (|T^) is D'^{x,Q'^) = Dg{x,Q'^)Dg{x,Q'^) 



Since the opacities used in (|13])-(0) for calculating D'^{x, Q^) were calculated using a Gaus- 
sian profile, we examined the sensitivity of D'^{x,Q^) to the choice of S{b±), so as to validate 
our normalization for the t = differential cross section which had been used for calculating 
the integrated cross section. The result of our calculations show that replacing a Gaussian 
profile with (||) leads to a change of which is less than 4% . We thus conclude that as far 
as screening corrections are concerned, we can safely rely on our previous estimations at t = 0. 

Eq. and Eq. can be viewed as the exchange of many "hard" Pomerons (gluon lad- 
ders). In the region of small \t\ we approximate this process by the exchange of one and two 
"hard" Pomerons (see Fig. The t-dependence of a single Pomeron exchange amplitude is 
proportional to ez^^^^*. The exchange of two "hard" Pomerons leads to an integration over 



where t = —q^. Hence, in Fig. ^ the first diagram is proportional to e^^^^^^ and the second 
diagram is proportional to e^^^^^^. 



We therefore approximate the differential cross section as the sum of these two exchanges 



I.e., 



= Al (e^^(^)* - .^ei^(^)f . (17) 

We fix our normalization factors A^^i and r] by the requirement that at t = 0, Eq. (|T2D and 
Eq. ([T7|) , once screening corrections as well as Fermi suppression are taken into account, match. 
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Comparing ( [I^ ) and (|T^ including these corrections we find: 



^(i.q") = l-r>(,x,e'), (18) 

where i^'/ is the Fermi motion normahzation factor. 

We calculate da jdt using (|T^ and (^) for low values of t, and compare the numerical cal- 
culations with the experimental data for elastic exclusive J/'?/' production [|1^. Note that our 



model, as presented in (plTD -fprsD, contains no additional parameters. Specifically, the normal- 
ization at t = is taken from the integrated cross section data and B is taken from section 
1^. Our results are shown in Fig. ^. The shaded areas in Fig. ^ correspond to the numerical 



calculations of (|T^ at each energy range given by Ref. [Q, and the dashed lines were calcu- 
lated by putting r/ = in It can be seen that a single-term exponential expression cannot 
describe the experimental data without modification to the overall normalization, which, as 
stated, is taken from rich data of integrated cross section. On the other hand, the reproduction 
of the data is reasonably good, when considering both terms of Eq. (P^). 

It is interesting to examine numerically the relative role of the first and second contributions 
of ([T7|). Squaring the ratio between the second term and the first term we have 7y2g2-^(^)l*l^ 
which is not only an increasing function of |t| but also an increasing function of W ^ since both 
B and ri increase with the energy. Hence, the second term becomes more and more important 
as |t| and/or W increase. According to our calculations, at t = the contribution of the second 
term to the differential cross section varies from 4% at = 50 GeV^ to 10% at = 150 GeVl 
At \t\ ~ Qo both the first and the second terms of ([T7|) are numerically larger than the value 
of da/dt, where the second term is 20(60)% from the first term at low (high) center of mass 
energy. 

As the theoretical uncertainties are difficult to estimate, we used the experimental errors to 
calculate the deviation of our simplified model from the data. The corresponding x^/n.d.f for 
the entire measured energy range is 0.7. As in [|l|, our parameters are strongly constrained by 
the experimental data. 



4 Differential Cross Section at Intermediate |t| 

The inelastic contribution to the cross section becomes more and more important as the mo- 
mentum transfer increases. The exponential t dependence (p!7|) was written using the elastic 
forward slope derived in section ^ which is compatible with the experimental data. In this 
section we would hke to suggest a simple expression for the inelastic differential cross section. 
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W =40 -54 GeV 



W =54 -66 GeV 




Figure 4: Elastic differential cross section at small t for different energy ranges. The shaded 
areas correspond to each energy range, as given by Ref. [^], the dashed lines 
correspond to a single-term exponential form, as further detailed in the text. 



10 



Figure 5: The proton vertex at —t > Qq. 



As stated in the introduction, hard processes can be calculated using the DGLAP equation, 
provided the scales of the interaction are such that the momenta on the gluon ladder are 
ordered. The production of a heavy vector meson can be viewed as the exchange of a gluon 
ladder between a nucleonic target and vector meson, which is a qq system of small transverse 
size (~ 1/My). For intermediate values of \t\, the momenta are still ordered on the ladder. 
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hence we argue that the DGLAP evolution continues to play a role up to |t| of the order of My 
Specifically, we can write the amplitude of the lowest rung of the ladder at the proton vertex 
(see Fig. 1^) using the two gluons proporgators A;f^ and /c^^. The corresponding amplitude is 
proportional to: 



2 



^2- 



(19) 



where ki = {q/2 + k), k2 = {q/2 — k) and the factor ki ■ k2 is implied by gauge invariance. 
From ([T9|) one can see that the dominant contribution to the evolution equation comes from 
large logarithms which appear when k"^ > = —t/4. For \t\ which is not too large, z.e.in the 
intermediate region, the DGLAP equation is still valid for calculating the production amplitude 
utilizing a gluon distribution function which evolves from —t/4 rather than from Qq. 



The physical interpretation of the above discussion is simple. At small \t\ (< l/rrip) the size 
of the target is roughly the proton size, the "evolution range" is large and the evolved gluon 
distribution consists of a large number of gluons. At higher values of momentum transfer the 
effective size of the target is of the order of l/|t|. This leads to a smaller number of emitted 
gluons, or in other words to a smaller numerical value of xG. 



From Eq. (|I2D we note that the differential cross section is determined by xG. The gluon 
distribution is related to the numerical value of the opacities, thereby to the screening correc- 
tions to the process. Thus, the main t dependence of the process comes from the scale from 
which the gluon distribution evolves. 

Based on the above discussion, we propose the following expression for the differential cross 
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W ^ 100 GeV 




-t [GeV2] 

Figure 6: Inelastic differential cross section in the intermediate region of Our calculations 
are presented, for two different values of the free parameter to, together with the 
experimental data of HI preliminary (squares), ZEUS 1995 (open triangles) and ZEUS 
preliminary (solid triangles) 



section: 

2 

^^i^^RA. = AlD\x,Q\ t) xG\x,q\ t) (20) 

where, 

xG{x,Q\t) = xgIx,-^-^], (21) 

and D(x,Q , t) is calculated by substituting ( pTD in (^ and carrying out the screening correc- 
tions procedure [Q] in both sectors. 

The justification for choosing the particular dependence of variables that the gluon distribu- 
tion depends on is to incorporate its depenence on t. Eq. (|2T|) is written in a typical dimensional 
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form, i.e.the second argument of the gluon distribution has dimensions of GeV^. This form 



is useful for most of the available parameterizations ||2^, |2^, ^ , where each parameterization 
evolves from a slightly different scale Qq. Since the DGLAP evolution is manifested by loga- 
rithms of Q'^/Qq, the gluon distribution function has a dimensionless form. Our ansatz may be 
better understood by considering a non-dimensional form of (0), namely xG{x, Q / {Ql+\t\/4:)) 
were it can be seen that as |t| increases, it dominates the scale of evolution whereas at small \t\ 
the gluon distribution function coincide with xG{x, Q/Ql). 



In (pOD, the coefficient is determined by a continuity requirement, that (|T^ and (pO| 



match at t = to, where to is our separation parameter. The results of the numerical calculations, 
for two different values of to, are shown in Fig. together with the HI and ZEUS data which 
have been read off the plot in |jl3| and [l^, respectively. As can be seen from Fig. ^ the large 



|t| preliminary data set of ZEUS [O] is below the HI preliminary data [113 . 



Our calculations of Eq. ([201) depends on the fitted value of our free separation parameter 
to- The upper curve shown in Fig. ^ corresponds to |to| = 0.54 GeV^, at which A^^ has been 
determined to reproduce the HI data with x^/n.d.f = 0.65. However, as the preliminary sets 
of HI and ZEUS data are consistently different, the task of finding a satisfactory description of 
both sets using the same parameters is not trivial. The best for the ZEUS data is obtained 
with |to| = 0.65 GeV^ (the lower curve of Fig. ^) where we calculated x^/n.d.f = 1.53 . Note that 
the difference between the two fitted values of to refiects the presumed difference between the 
HI and ZEUS sets. This is not necessarily a genuine difference and may just be a consequence 
of different cuts applied at the diffractive proton vertex, by the two collaborations. 



At first sight, it looks as if ( poD may affect the results of our previous publication ||l]|, inter 
alia on account of the t-dependence of screening corrections at high |t|. However, this effect 
is compensated by the value of the forward slope, which we calculate using the dipole electro- 
magnetic form factor, as opposed to an exponential form factor used in Indeed, we found 
that we still reproduce all the available experimental data. The compensation is so accurate, 
that the two curves practically overlap, and we have the same values as given in 



5 Summary and Conclusions 

In this paper we have explored the momentum transfer dependence of exclusive production of 
J/ip vector mesons. By using a Fourier transform of an electromagnetic form factor as the 
profile function in the impact parameter space, we calculated the forward differential cross 
section slope. 

We used the expression derived for the slope and calculated the differential cross section for 
the elastic production of J/ ip by using a simple exponential expression approximating one and 
two Pomeron exchange. 
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For the inelastic process, we calculated the differential cross section as a function of t, based 
on a somewhat naive picture of the interaction. In which we suggested replacing the argument 
of the gluon distribution so as to obtain a decreasing function of which coincides with pQCD 
calculations at small \t\. 

Our conclusions are: 

1. The value of the B slope is sensitive to the choice of the profile function S{bj_). Our 
calculations of B reproduce the experimental data and are not sensitive to the choice of 
the PDF. 

2. The fact that the screening correction formalism with a dipole-type b± dependence {FfUpoie) 
reproduces both the value and the energy dependence of the t-slope, is important for 
the understanding of the energy dependence of the "hard" Pomeron trajectory. The 
contributions of the screening corrections increase with energy and therefore ctg^j, the 
effective slope of the "hard" Pomeron trajectory, should also increase with energy. This 
result is particular to a screening corrections approach and is not obtained in non-screened 
pQCD calculations. 

3. The screening corrections, which we express through the damping factor D are not effected 
by a change in the h dependence of the opacities. 

4. An exponential form of the differential cross section, which is in agreement with the 
experimental fits to the forward slope, is suitable for reasonable description of the elastic, 
low experimental data. This form consists of the first two terms of a multi Pomeron 
exchange, which are supplemented by screening corrections, as well as corrections due to 
Fermi motion of the quarks within the charmonium system. This approximation is also 
valid for describing the integrated cross section. 

5. At intermediate values of inelastic data are well reproduced by pQCD which is based 
on DGLAP evolution. The value of \t\ at which the inelastic processes are comparable 
with elastic processes is —to ^ 0.5 GeV^. 

6. At intermediate values oi\t\ the screening corrections dependence on t reflects the decrease 
of the gluon distribution function. 
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